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Abstract 

Applying domain decomposition to the lattice Dirac operator and the associated 
quark propagator, we arrive at expressions which, with the proper insertion of ran- 
dom sources therein, can provide improvement to the estimation of the propagator. 
Schemes are presented for both open and closed (or loop) propagators. In the end, 
our technique for improving open contributions is similar to the "maximal variance 
reduction" approach of Michael and Peisa, but contains the advantage, especially 
for improved actions, of dealing directly with the Dirac operator. Using these im- 
proved open propagators for the Chirally Improved operator, we present preliminary 
results for the static-light meson spectrum. The improvement of closed propagators 
is modest: on some configurations there are signs of significant noise reduction of 
disconnected correlators; on others, the improvement amounts to a smoothening of 
the same correlators. 
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1 Introduction 



In recent years, there has been a number of techniques developed for improving 
estimates of quark propagators on the lattice (see, e.g., Refs. [1,2,3,4]). The 
main idea is to try to devise an estimator formulation which allows one to use 
all (or at least many) lattice sites as source locations for the quarks, rather 
than just having a fixed location, which is traditionally the case. If successful, 
one then has an estimate of the quark propagators from anywhere to anywhere 
in the lattice (typically called "all-to-all" propagators). 

In the following sections we add our voices to the crowd and present our 
own method, which is based, foremostly, upon domain decomposition, or more 
specifically, the decomposition of block matrices. We end up with two different 
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types of estimators: one for "open" propagators between two domains of the 
lattice (or "half-to-half" ) and one for "closed" propagators within one of the 
domains. 

We first show how we devised our method through the consideration of random 

sources. Later we present first results for the two types of estimators and some 
of their possible applications: static-light mesons for the open propagators and 
disconnected correlators for the closed. 



2 The method 



Using a set of random sources, XjfBb^ placed at all sites of the lattice, one 
can determine a corresponding set of solution vectors (we use a summation 
convention for repeated indices throughout this paper): 

vZ^a-MrljPbXm in^l,...,N). (1) 

Since we know that lim^v^oo jfXj/SbXklc — ^jk^P'y^bci we can construct an esti- 
mate of the full propagator: 

iaak-yc ~ ji^ liaaA.k'yc ' \ ) 

This is, however, a rather noisy estimator and we map out some improvement 
schemes in the following sections. 

The main feature of our approach involves the consideration of independent 

regions of the lattice, or domain decomposition, a technique used previously for 
"maximal variance reduction" (MVR) of pseudofermion estimators [1]. Similar 
considerations may also be used for improving dynamical lattice updates [5]. 

We can think of the lattice as a disjoint union of two regions. The full Dirac 
matrix can then be written in terms of submatrices 

_ / Mn Mi2 
\^ M21 M22 

where Mn and M22 connect sites within a region and M12 and M21 connect 
sites from the different regions. (For simplicity, we suppress color, spin, and 
site indices in the following.) Regardless of the shape or nature ^ of the regions, 
a similarity transformation is all that is needed to reach this form. We can 

^ For example, different regions in color or Dirac space. 
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also write the propagator in this form: 



M-^ = P = 




(4) 



2.1 Open contributions 



It is helpful to imagine reconstructing the sources in one region, Xi, from the 
solution vectors everywhere, and to separate the contributions from the 
different regions: 

X^ = M?7" = Mn77r + Mi2772"- (5) 
If we now apply the inverse of the matrix within one region, we have 



(6) 



This can be solved for 77" and substituted into the original expression for the 
noisy estimator of the propagator between the two regions: 



12 



-n.2 ~ -^^1 X2 



1 

TV 



1 {m^:m,,7i^) , 



N 



(7) 



where in the last line we eliminate the first term due to the fact that we expect 
lim^v^oo XiX2^ = 0- This is a crucial step, for here we cut out of the calculation 
what amounts to being only noise. It does not come for free, however, since 
we must perform the additional inversion within the subvolume 1. Writing out 
the full expression, we obtain 



Pi2~-^MniMi2Px2xf 

= -Mf/ M12P22, (8) 

where the second line is an exact expression, showing that one can relate 
elements of different regions of P = via the inverse of a submatrix of M. 
(We do not pretend to have derived something new here; after all, P22 is the 
Schur complement of Mii. We only wish to emphasize the useful connection 
with random-source techniques.) Again, the lesson learned up to this point is 
that we need no sources in one of the two regions. The story does not end 
here, however. 
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Looking again at Eq. (8), one can see that we need not make the approxima- 
tion P22 ~ jfPX2X2^ ■ Instead, we can place the approximate Kronecker delta 
between the M12 and P22: 

Pl2~-^Mr/Mi2X2X2^P 

~-^(Mr/Mi2X^) (75^75%^)^ 

^-^V'r</'2^ (9) 

where we have used the 75-hermiticity of the propagator (this is done only for 
convenience since we could just as well work with P^X2 '^2 )- Fiora the 

next to last line, one can see that the signal only rises from terms where the two 
components of the X2 s the same. However, unlike the naive estimator, Eq. 
(2), where there is only 1 such term giving a signal-to- noise of ~ l/F^/^, here 
there are many: for source points, the signal-to-noise is ~ V^/ {V^ — Vy^Y^'^ ~ 
1. Terms where the components of the sources are not the same can still be 
eliminated by "dilution" of the original source vectors, that go into the 
ipi, the 02 ) for greater noise reduction, both simultaneously. Probably more 
important than these considerations, however, is the fact that, for most of the 
propagators between the regions, the random sources are kept far from the 
end points. Also, one can use all points in one region for the source and all 
points in the other region for the sink. 

The ideal domain decomposition for quark propagators which contribute to 
connected diagrams is then to use two unequal volumes, one containing a few 
more time slices (those of the sources x) than the other. Ideally, the centers of 
the two sets of source time slices should be separated by T/2. The number of 
source time slices is dictated by the lattice Dirac operator since the x's should 
be placed on all time slices which communicate with the other region via one 
application of M. For Wilson and Fixed-Point (FP) [7] quarks, this is just 2 
time slices, 1 adjacent to each boundary. For Chirally Improved (CI) quarks 
[6], which we use here, 4 are necessary (see Fig. 1). For the Asqtad action [8], 
6 are needed due to the presence of the Naik term. While for Overlap fermions 
[9] , it might be best to use equal volumes for the two regions since the sources 
will have to cover one region entirely. But we stress that for all of the above, 
one is free to dilute further: e.g., by inverting the sources on the different time 
slices separately. 

For our first attempt of using this method, we do not choose the ideal de- 
composition. We use equal volumes for the two regions and place sources next 
to the boundaries in both regions (see Fig. 2). Although this choice may not 
be ideal, we do perform inversions for each spin component separately (spin 
dilution). Note also that, since our sources occupy all relevant time slices sur- 
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Fig. 1. Depiction of the ideal decomposition and the sources (x) which are needed 
to construct the estimate of the CI quark propagator between regions. 
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Fig. 2. Depiction of one of the boundaries and the surrounding sources (x) which 
we use to construct the two estimates of the CI quark propagator between regions 
of equal volume. Colors and shapes indicate which source region contributes to the 
signal in the resultant vectors. The ip^s are calculated using only one of the source 
regions, while the 0's use both. 

rounding the boundaries, we can actually obtain two independent estimates 
of the quark propagator between the two regions: 



Pi 



1 



12 



(10) 
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So our method is very similar to that of MVR, except for the fact that we 
can work directly with M, instead of Af"'"M. This gives the current method 
three advantages: First, it is less problematic to invert M since it has a bet- 
ter condition number than M^M [1]. Second, implementing the method is 
straightforward. One needs only to limit the range of application of M when 
performing inversions in the subvolumes and when multiplying by the matrix 
between regions. Otherwise, existing routines remain imchanged. Third, the 
sources need only occupy enough time slices to connect them with the other 
region via M, rather than M^M. These are the same number of sources for 
Wilson-like operators, where M^M, like M, only extends one time slice. How- 
ever, for many other improved operators (like CI and FP) this can reduce the 
number of necessary source time slices by a factor of 2. 



2.2 Closed contributions 



Working with the propagator within one of the regions, say 1, and following 
a similar derivation as in the previous section, one obtains an expression po- 
tentially useful for estimating quark propagators which return to the same 
region: 

=Mr/-Mr/Mi2P2i. (11) 

So once again, through the consideration of random sources, we arrive at an 
exact expression (and again, one which is nothing new). Now, combining the 
expressions for Pn and P21 (= TsPiWs); we arrive at the relation: 

Pii = Mfi^ + MniMi275 (Mf/Mi2P22)%5 ■ (12) 
Inserting our random sources into this expression gives 



Pn - Mr/ + ^ (Mr/Mi275X2) {l.M^i^ M^^Px' 



N 

^M-V^CiX^ (13) 

Note that we put no region index on the second %, indicating that for this 
resultant vector (^1) we wish to use sources initially placed everywhere on the 
lattice. The advantage of expression (13) may not be immediately clear since 
it still contains the explicit appearance of M^^. 

This seeming hindrance can be remedied by considering a very small volume 
for region 1. Performing this "highly reduced" inversion exactly, we hope to 
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Fig. 3. Depiction of the boundary and the surrounding region of sources which 
are needed to construct the estimated CI quark propagator to and from the central 
point. In this figure, region 1 is just one site (we also consider a larger region 
which includes nearest neighbors) . The red squares indicate which region of sources 
contributes to the signal in the resultant vectors (Ci, ^i) of Eq. (13). 

find a significant gain in the signal-to-noise ratio: The first term in the above 
expressions will be exact and the second term, as compared to its naive esti- 
mate, may be improved by a factor of as much as ~ '^x^'^^ where Vy^ is now 
the volume of sources in region 2 which connect to region 1 via M. So if the 
volume of region 1 is kept small enough and the lattice Dirac operator con- 
nects each site to many others, there may be an advantage to calculating the 
M{i^s exactly, as opposed to simply inverting more sources. 

Here, we propose to use spin-diluted, but full-volume sources for x, so that we 
can, in the end, use all sites as the starting and ending point of the propagator 
Pxx with only original full- volume inversions: the Px" iii ^i- We can consider 
the smallest, symmetrical volume for region 1, the point itself {Vi = 1; see Fig. 
3), in order to reduce the amount of work needed to calculate M^i, which we 
need for each point in the lattice. With this choice, at most one needs to invert 
^ 12 X 12 matrices. Since = 128 for the CI operator in this scenario, we 
hope that this small amount of extra work may be worth the effort. In the 
end, however, we actually use a larger region 1, including nearest neighbors 
(Vi = 9), but we avoid some extraneous work by only inverting Mu upon 12 
sources located at the central point. 



3 First results 

In this section, we present some first results obtained with our estimated quark 
propagators. We first use the open propagators for the light quark in heavy- 
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Table 1 

Static-light meson operators. 



oper. 




QO(r,A)g' 


S 


o-,i- 




P_ 


0+,l+ 




P+ 


l+,2+ 


Q (71^1 - I2D2) q' 


D± 


l-,2-,3- 





light meson correlators, where we use the static approximation for the heavy 
quark. We also show the effects of using our closed propagators in disconnected 
correlators. 



3.1 Static-light mesons 



For our meson source and sink operators, we use bilinears of the form: 

QO{V,D,,D\S)q, (14) 

where 5 is a gauge-covariant (Jacobi) smearing operator [10] and D is the 
covariant derivative. For our basis of light-quark spatial wavefunctions, we use 
three different amounts of smearing and apply 0, 1, and 2 covariant Laplacians 
to these: 

q'= S^q, D^Suq, D^D''S,,q, (15) 

where the subscript on the smearing operator denotes the number of smearing 
steps; all are applied with the same weighting factor of k,,5„j = 0.2. So wc have 
a relatively narrow, approximately Gaussian distribution, along with wider 
versions which exhibit one and two radial nodes, due to the application of 
the Laplacians. We point out that, thus far, we have not altered the quantum 
numbers of the meson source since both the smearing and Laplacians treat all 
spatial directions the same (i.e., they are scalar operations). In order to create 
mesons of different quantum numbers, we use these light-quark distributions 
together with the operators shown in Table 1 (see, e.g., Ref. [1]). 

Inserting the estimated and static propagators in the meson correlators we 
have 



Qj{t) = (0 {QO,q)t {qO,Q)o 



ETr 



1 + 74 



t-i 



l[ul{x + il)0,P,^,i^,0. 



i=0 



(16) 



{U} 



The static quark is propagated through products of links in the time direction 
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and has a fixed spin (1 + 74)/2. The estimated propagator P^j^ia x the 
form of Eq. (10). Thus, all points within region 1 of them) can act 

as the source location x, just so long as t is large enough to have the sink 
location x + tA in region 2. Note that we now have subscripts on the source 
and sink operators to denote which light-quark distribution is being used. We 
create all such combinations and thus have a 3 x 3 matrix of correlators for 
each of the operators in Table 1. 

Following the work of Michael [11] and Liischer and Wolff [12], we use this 
cross-correlator matrix in a variational approach to separate the different mass 
eigenstates. We must therefore solve the generalized eigenvalue problem 

C(t)iT(^)=A«(t,to)C(to)i^'\ (17) 

in order to obtain the following eigenvalues: 

X^''\t,to) oc e-*^^[l + 0(e-*^^^)] , (18) 

where is the mass of the A;th state and AM^ is the mass-difference to the 
next state. For large enough values of t, each eigenvalue should then correspond 
to a unique mass state, requiring only a single-exponential fit. 

Before solving the above eigenvalue problem, we check that our cross-correlator 
matrix is real and symmetric (within errors). We then make it symmetric 
before solving Eq. (17) and look for regions of t where the eigenvalues exhibit 
single masses and the corresponding eigenvectors, v^''\ are stable (each of these 
provides what we call the "fingerprint" of a state and if it holds steady, we 
have better reason to believe that we are looking at a single state). The results 
should be independent of the normalization time slice to and we make sure 
that this is so. For our final results, however, we use the value oi to/a — 1. 

This variational approach has seen much use recently in lattice QCD, espe- 
cially for extracting excited hadron masses and we point the reader to the 
relevant literature in [13]. 

We create our cross-correlator matrices on two sets of gluonic configurations: 
100 quenched and 74 dynamical, each with 12^ x24 lattices sites. The quenched 
configurations have a lattice spacing oi a ^ 0.15 fm (a~^ f« 1330 MeV) and 
a spatial extent of L f« 1.8 fm. The dynamical set [14] has 2 flavors of CI 
sea quarks (with M^^sea ~ 500 MeV), a 0.115 fm (a'^ ^ 1710 MeV), and 
L 1.4 fm. We use 12 random spin-color vectors as sources for the light-quark 
propagator estimation. Spin-diluted, this gives us 48 separate sources for the 
inversions (one in the full volume, 0, and two in the subregions, ip; see Eqs. 
(9) and (10)). We perform inversions for 4 different quark masses: aruq — 0.02, 
0.04, 0.08, 0.10. 

After extracting the eigenvalues, we check for single mass states by creating 
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Fig. 4. Effective masses for the static-light mesons on the quenched configurations. 
aniq = 0.08, ~ 1330 MeV, L ~ 1.8 fm. The horizontal lines represent M ± 
fit values for the corresponding time ranges. Dashed lines indicate fits for which we 
adjust the minimum time for systematic error estimates. 



effective masses: 



aM 



(fc) 



eff 



In 



/ AW(t) 



(19) 



A representative sample of these, along with single-elimination jackknife er- 
rors, are plotted against time in Figures 4 and 5. In each case one finds values 
from the first two or three eigenvalues. The horizontal lines signify the M±(Tm 
values which result from correlated fits over the corresponding range in time 
(dashed lines denote fits where we also later increased the minimum time of 
the fit in order to estimate systematic errors). We require that at least three 
effective mass points display a plateau (within errors) and that the eigenvec- 
tors remain constant over the same range before we perform said fits. Figure 
6 shows the eigenvector components for the quenched first-excited P-wave 
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Fig. 5. Effective masses for the static-light mesons on the dynamical (Mjr^sea ~ 500 
MeV) configurations, aruq = 0.08, a-^ ^ 1710 MeV, L ^ 1.4 fm. The horizontal 
lines represent M ± ctm fit values for the corresponding time ranges. Dashed lines 
indicate fits for which we adjust the minimum time for systematic error estimates. 

(2P+), a state displaying a rather jumpy effective-mass plateau. There is some 
slight variation in the central values over time, but when considering the fine 
scale of the plots and the errors, we are confident that we are dealing with a 
single state. 

Performing fits for all quark masses, we next take a look at the mass splittings 
(M — Mis) as a function of the quark mass. These are plotted in Figure 7, 
along with the chirally extrapolated results [rriq — ^ 0). We use simple linear fits 
to perform these extrapolations (as well as for the interpolations to rrig = iris) 
since there appears to be little scaling in these plots. Both sets of results display 
a clear ordering of alternating parities with increasing mass. The quenched 25" 
and 1D± states are close, with the latter being slightly higher in mass. This is 
expected since the D± represents an average of higher spins than the S and the 
heavy-quark spin interactions which would "mix" the purely orbitally excited 
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Fig. 6. The three eigenvector components for the quenched 2P+ state over the time 
range of the fit. 
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Fig. 7. Physical mass sphttings (M — M15) as a function of the quark mass for 
the quenched (left) and dynamical (right) lattices. The vertical lines denote the 
chiral limit (m^ — > 0) and the strange quark mass (m^). The left-most values are 
the (linear) chiral extrapolations. All errors are only statistical for the chosen fit 
ranges. 

D-waves with the ground-state S-wave are absent in the static approximation. 
Overall, the results for the dynamical configurations are poorer. This is no big 
surprise since here we have not only fewer configurations, but also a smaller 
physical volume (L 1.4 fm). This fact also makes it difficult, at least at the 
present stage, to discern quenching and finite- volume effects. However, the one 
marked difference, the jump in the 35* mass, is quite likely due to the smaller 
volume. 
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Table 2 

Mass splittings for our B mesons {niq —i- 0). The first number in parentheses is 
the statistical error. The second set (if present) are the additional systematic errors 
which result from adjustments to the minimum time of the fit. 



state M- Mis (MeV) 

Nf = 0,L^ 1.8 fm Nf = 2,L^ 1.4 fm PDG [15] 



2S 
3S 


0", 

0", 


1- 
1" 


712(14) 
1265(40)(i;3o) 


717(69) 
1640(44)(lio) 




1P_ 
2P_ 


0+, 
0+, 


1+ 
1+ 


350(46) 
971(49)(1^0) 


243(51) 


384(8) 


1^+ 

2P+ 


1+, 

1+. 


2+ 

2+ 


446(15) 

1028(28)(;j[]") 


341(82) 

1001(80)(;^^") 


384(8) 


1D± 

2D± 


l-,2 

l-,2 


-,3- 
-,3- 


808(27)(l°o) 

1183(97)(tl^0) 







In Tables 2 and 3 we present the results for the chirally extrapolated {B 
mesons) and strange-quark-mass interpolated {Bs mesons) mass splittings. 
We include the statistical errors from the fits in the first set of parentheses. 
For fits where the effective-mass plateaus are not immediately clear (e.g., fits 
represented with dashed lines in Figures 4 and 5), we move the minimum 
time of the fit range out by one to two time slices and observe the subsequent 
changes in M ± ctm, as compared to the previous values. The differences from 
the old values are reported as systematic errors; these appear in the second 
set of parentheses. 

It is interesting to note that, in moving from the quenched lattices to the 
smaller, dynamical lattices, the M2S — Mis mass sphtting remains unchanged. 
Unless there is an odd counter-balancing of finite- volume and quenching effects 
between the lattices, we may conchidc that our values for M25 — Mis are 
reliable. The same is true for the M2P — Mis splitting, which is a bit strange 
since Mip — Mis already shows a difference, its value becoming especially low 
for the IP-, when compared with the result from experiment (the quantum 
numbers, J^, for the experimentally observed, excited B and Bs mesons have 
yet to be confirmed, so we place them in both the 1P__ and 1P+ rows). For our 
Bg mesons, all Mip — Mis splittings appear to be too small. Better statistics 
and larger volumes are needed if we are to clearly resolve these matters. 

One thing is clear though: due to the improvement of the light-quark propaga- 
tor estimation, and our subsequent ability to use half the points of the lattice 
as source locations, we have greatly improved our chances of isolating excited 
heavy-light states. In an earlier study [16] of heavy-light mesons using wall 
sources on the same quenched configurations, we were barely able to see the 
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Table 3 

Mass splittings for our Bg mesons (m^ = m^). The first number in parentheses is 
the statistical error. The second set (if present) are the additional systematic errors 



which result from adjustments to the minimum time of the fit. 


state 




M 


- Mis (MeV) 








Nf = 0,L^ 1.8 fm 


Nf = 2,L^ 1.4 fm 


PDG [15] 


2S 


o-,i- 


675(10) 


665(45) 


- 


3S 


o-,i- 


1220(30)(;20) 


1560(45)(l?^o) 


- 


1P_ 


0+,l+ 


384(20) 


330(34) 


448(16) 


2P_ 


0+,l+ 


923(30)(li0) 






1^+ 


l+,2+ 


424(10) 


363(55) 


448(16) 


2P+ 


1+.2+ 


993(20)(;^^") 


930(75)(i;|,) 




1D± 


l-,2-,3- 


773(17)(l°o) 






2D± 


l-,2-,3- 


1188(68)(t^o'°) 







25* state, let alone the excited states in any other operator channel. Also, there 
we used NRQCD for the heavy quark; this should only boost the signals since 
the heavy quark can then "explore" more of the lattice through its kinetic 
term. It is obvious, however, that we have much better signals now since we 
are able to see excited states in every channel (25', 3S, 2P_, 2P+, and 2D±) 
on the quenched lattice. 

3.2 Disconnected correlators 

As a preliminary testing ground for our closed propagators, we take a look at 
some first results for the disconnected contributions to pseudoscalar (J^ — 0~) 
meson correlators: 

D{t) = Tr(75Px,to;^,to)Tr(75^y-to+t;y,to+t) • (20) 

to Ay 

Again, we use 12 random spin-color sources (initially placed everywhere on the 
lattice), spin-dilute them, and perform inversions (Px") a quark mass of 
arUq = 0.02 on the quenched configurations. We then condition these "naive" 
estimates via Eq. (13) using the central point and its nearest neighbors (NN) 
as region 1 (the calculation of all the Mj^^'s for this sized region on a single 
configuration takes less than a day on a PC). 

In Figure 8 we compare results obtained via the naive and improved estimators 
on two different quenched configurations. The errors are estimated via the 
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Naive and improved disconnected pseudoscalar correlators (with 



aruq = 0.02) on two different quenched configurations. Errors result from a sin- 
gle-elimination jackknife procedure over the 12 random spin-color sources. 

single-elimination jackknife subsets of the 12 random sources. Looking at the 
result for the configuration on the left, we can see significant reduction of the 
errors over many time separations. This is not the case for all configurations, 
however, as one can see on the right, where the errors are comparable, if not 
slightly larger for the improved version. For both cases shown here (and in 
fact for all configurations) the central values for the improved method follow a 
smoother curve. This should be no great surprise since the improved estimator 
uses sources on neighboring time slices (see Figure 3), whereas the naive one 
does not. So the improved version should smoothen out some of the remaining 
fluctuations over different t values. 



An important thing to note here is that, despite any of the improvement 
which we may gain from the smaller errors on some configurations and the 
smoothening of the curves, the error for the ensemble average of D{t) will 
be dominated by the limited number of gauge configurations (i.e., this is a 
"gauge-limited" quantity). One can see from the figure that the fluctuation 
which comes from switching configurations is as big as, or bigger than, the 
jackknife errors from the sources. For some other configurations, the jump in 
the corresponding D{t) values is much larger. These are perhaps configurations 
with large values of topological charge, Q; after all, the integrated disconnected 
pseudoscalar correlator is related to the square of this quantity [17]: 

T.m^—,, (21) 

t 

where the relation is only approximate here since we use only chirally improved 
quarks. 

Although we see no great overall improvement for these light, disconnected 
pseudoscalars, this is a tough testing ground. It remains to be seen how well 
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these improved closed propagators will perform for other quantities where only 
averages of single loops over single time slices are needed: e.g., ss loops within 
hadron correlators to measure strangeness content, an application which we 
plan to look into in the near future. We also point the interested reader to 
other improvement schemes for closed contributions: see, e.g., [18,19]. 



4 Summciry 

We have presented a method for improving the estimation of quark prop- 
agators in lattice QCD. The improvement we obtain for open propagators 
between equal subvolumcs of the lattice is significant. The boost in statistics 
from being able to use up to half of the lattice sites as starting locations for 
our static-light meson correlators allows us to extract a number of radially 
and orbitally excited states {2S, 3S, IP, 2P, ID, 2D). 

As we have pointed out, our method is similar to the maximal variance reduc- 
tion approach [1] . However, since we can work with the Dirac matrix M, rather 
than M^M, our method may be better suited for highly improved actions. 

Although the first results for our improved closed propagators arc not very 
promising, it appears that we have chosen a very difficult application (light, 
disconnected, pseudoscalar correlators) for which to test them. We reserve 
final judgment on their usefulness for the future, when we measure their per- 
formance in other physical quantities. 
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5 Appendix A 

In Section 2, we derived the exact relations (8) and (11) for the blocks of the 
quark propagator. In those derivations we relied upon the following property 
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of the random vectors , n — 1, . . . , N: 



^2L{h>^?)-''r (22) 

But there also exist other possibihties to derive the same results. 

As an alternative, we show here how one can arrive at Eq. (8) using the path 
integral over the fermionic degrees of freedom. 

The Dirac propagator, P — M"^, can be represented by a path integral in the 
following way 



Pij = I / [Di^] [Di,]^Pii,^e-'f , (23) 



where 



Z^J[Dij][Dij]e-'f, (24) 
is the partition function, Sf is the fermion action 

Sf^jj^Mi/;, (25) 
and [Dip] = H/t dipk is the integration measure. 

The first step in our alternative approach is to split the fermion action into a 
sum of four separate contributions: 

Sf = IplMulpi + V^2M22V'2 + i^lMi2lp2 + i^2M2l'lpl. (26) 

The leading two contributions are only connecting the fermion fields within 
the regions 1 {Ri) and 2 (-R2), while the other two contributions provide the 
connections between the regions. In a similar way, one can also divide the 
integration measures in Eqs. (23) and (24). 

When one now considers the special case of Eq. (23) when i e i?i and j G R2 
and uses expression (26) the propagator becomes 
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Pi2 = ^ J[Di;][Dij]^,ij2e-''f 



Z 

det{Mu] 



[D^2][DiJ2]M-M^,'MuiJ2)e-^'''''^' 



(27) 



where Mu and M22 := M22 -MaiMfi^Mia are the diagonal blocks of a decom- 
posed version of M with O's as off-diagonal blocks. Therefore, Z = det{M) = 
det{Mii)det{M22) and the ratio becomes 1. For the next to last step, we com- 
pleted the square in the exponent and performed a shift of the integration 
measure. 

Since M22 is the Schur complement of Mu, one can show that in fact P22 = 
M22^. This leads us to our final result 
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